arXiv:1504.02156vl [astro-ph.HE] 8 Apr 2015 


Can black-hole neutrino-cooled disks power short gamma-ray 

bursts? 

Tong Liu 1,2,3 ’ 1 2 3 4 ’*, Yi-Qing Lin 5 6 , Shu-Jin Hou 6,2 , Wei-Min Gu 1,4 
*tongliu@xmu.edu. cn 

ABSTRACT 

Stellar-mass black holes (BHs) surrounded by neutrino-dominated accretion 
flows (NDAFs) are the plausible candidates to power gamma-ray bursts (GRBs) 
via neutrinos emission and their annihilation. The progenitors of short-duration 
GRBs (SGRBs) are generally considered to be compact binaries mergers. Ac¬ 
cording to the simulation results, the disk mass of the NDAF has been limited 
after merger events. We can estimate such disk mass by using the current SGRB 
observational data and fireball model. The results show that the disk mass of a 
certain SGRB mainly depends on its output energy, jet opening angle, and central 
BH characteristics. Even for the extreme BH parameters, some SGRBs require 
massive disks, which approach or exceed the limits in simulations. We suggest 
that there may exist alternative magnetohydrodynamic processes or some mech¬ 
anisms increasing the neutrino emission to produce SGRBs with the reasonable 
BH parameters and disk mass. 

Subject headings: accretion, accretion disks - black hole physics - gamma-ray 
burst: general - neutrinos 


1 Department of Astronomy and Institute of Theoretical Physics and Astrophysics, Xiamen University, 
Xiamen, Fujian 361005, China 

2 Key Laboratory for the Structure and Evolution of Celestial Objects, Chinese Academy of Sciences, 
Kunming, Yunnan 650011, China 

3 State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of 
Sciences, Beijing 100190, China 

4 SHAO-XMU Joint Center for Astrophysics, Xiamen University, Xiamen, Fujian 361005, China 

5 School of Opto-electronic and Communication Engineering, Xiamen University of Technology, Xiamen, 
Fujian 361024, China 

6 College of Physics and Electronic Engineering, Nanyang Normal University, Nanyang, Henan 473061, 
China 



2 


Introduction 


Gamma-ray bursts (GRBs) are the most powerful electromagnetic events in the uni¬ 
verse, which are sorted into two catego ries, i.e., short- and long -duration GRBs (SGRBs 


and LGRBs, see iKouveliotou et al.l Il993t) or type I and II GRBs (IZhandl2006t IZhang et al. 


2007aj). Their progenitors are considered to be mergers of two compact obje cts, i.e. , tw o 


neutron star s (NSs) or a black hole (BH) and a NS (for reviews, see, e .g., iNakarl 120071: 


Berger! 1201411 . and collapses of massive stars (e.g., IWooslev fe Bloom! 120061 for reviews), re 


spectively. For the interpretations on the gamma-ray and afterglow emission of G RBs, the 


widely accepted. The popular models on the central engines of GRBs are either a rotating 

stellar BH surrounded 

)V a hvoeraccretion disk (e.g., Paczvnski 1991; Naravan et al. 

1992; 

MacFadven & Woosley 

1999h or a cuickly rotating magnetar (or protomagnetar, e.g., 

Usov 


1992 


Metzger et al. 


2011 


Lii et al. 


2015). 


Two mechanisms have been proposed to power GRBs if a hyperaccretion disk exists in 
the center of GRBs, i.e., neutrino emission and annihilation, and_magnetohydrodynamic pro¬ 


cesses, such as Blandfor d-Znajek (BZ) mecha nism flBlandford fe Znaiekl 1197711 and episodic 


magnetic reconnection (Yuan fe Zhan g 2012). For the former mechanism, neutrino anni¬ 


hilation can produce a relativistic electron-positron outflow, which is considered as the 
progenitor of the fireball to power a GRB. The most probable model to launch a large 
number of neutrinos is a geometrically and optically thick neutrino-cooled hyperaccretion 
disks, named as neutrino-dominated accretion flow (NDAF), whose typical characteristics 
are extremely high accretion rate and neutrino-cooling process. In the inner region of 
the NDAF, the main components are the electrons, free neutrons and protons, the den¬ 
sity and temperature are very high (p ~ 10 10 — 10 13 g cm” 3 and T ~ 10 10 — 10 11 K), 


ir -o <7 — JT — —-— j | 

neutrino and antineutrino radiation (see, e.g., Popham et al. 

1991 

): 

3i Matteo et al. 

2002 

Kohri & Mineshige 2002; 

Rosswog & Ramirez-Ruiz 

20 

02; Kohri et al. 

2005; Lee et al 

2005 

Gu et al. 2006); Chen & Belo 

mrodovj 200?! 

Liu et al. 

200A 

kawanaka & Mineshige 

2C 

07 

Zalamea & Beloborodov 

2011 

Janiuk et al. 

20131: I 

kawanaka et al. 

20131 Xue et al. 2013 

)• 


Two factors should be considered in the calculation of the neutrino luminosity and 
annihilation luminosity, which are the structure and components of the NDAF and the de¬ 
scription of the relativistic neutrino propagation. Xue et ah (2013) investigated the global 
solutions of the radial structure and components of the NDAF in the Kerr space-time of the 
BH with detailed neutrino physics and nucleosynthesis processes. The results show that the 
gas pressure and the neutrino cooling are always dominant in the inner region for the high 
mass accretion rate, and the major components of the inner, middle, and outer regions are 
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the free nucleons, 4 He, and 56 Fe, respectively. Importantly, they noticed that the radiative 
neutrinos mainly come from the inner region of the disk, and the neutrino emission rate less 
depends on the description of the microphysics, as well as o t her studies of NDAF m odel (e.g., 
Popham et al. 1999 : Di Matteo et al. 2002 : Liu et al. 2007 : Kawanaka et al. 2013 1. Even for 
the d iscussions on the vertical structure of th e NDAF, the similar solutions are presented 
(e.g., Liu et al. 2008 . 2010 . 2012a . 2013 . 2014). Thus the main problem is how to pre cisely 
calculate the neutrino an nihilation processes. iBirkl et al.l (1200711 . iKovacs et al.l (l2011al) . and 
Kovacs fe Harkol (120 11b) analyzed the influence of general relativistic effects on the neutrino 
annihilation efficiency, which has a prominent inc rease co m pared with the New tonian calcu¬ 
lations. Based on the geodesic-tracing method, Zalamea & Beloborodov (20111 also studied 
annihilations via tracing the neutrino track. 


For SGRBs, Eichler et a 


(1989) proposed that the mergers of two NSs might be the 


candidates. Ruffert fc .Tanka f 1998b simulated three-dimensional Newtonian hydrodynamical 
solutions of the merger events of two NSs with mass ~ 1.6 M & . There might survive a 
disk ~ 0.1 — 0.2 Mp) surrounding a BH ~ 2.5 M & . Furthermore, Paczvhsk i (|199ll) and 


Nara van et al. (Il992f l presented that the merger of a NS and a stellar-mass BH can also 
produce a SG RB. In simulations, the fragments of the NS can form a more massive disk, 


(L5Mq (e.g., iKluzniak &; Ledll998t iLee &; Kluzniaklll999l: IPopham et al.lIl999c iLiu et al. 
2012b|). In the past several years, the massi ve NSs, i.e., ~ 2 have been discovered i n 
binaries, which accompany the white dwarfs (IDemorest et al.l l2010l: lAntoniadis et al.l 120131) . 
Yet we cannot neglect the possibilities that the massive NSs exist in the BH-NS or two NSs 
binaries. However, the mass of the disk is still much smaller than 1 M & with the logical 
conjecture. So it begs a question: can annihilati ons o f neutrinos from NDAFs owning such 
disk masses power all the observed SGRBs? Fan &l Wei (20111) investigated the disk mass in 
the center of SGRBs with the fixed values of the BH mass and spin. They found that nearly 
half SGRBs are suitable for the results of the above simulations. We further consider that 
the annihilation description, intact samples with prompt emission and afterglow properties 
of SGRBs, and reasonable ranges of the BH parameters should be fully included to answer 
the above question. 


In Section 2, we describe the physical processes from the neutrino annihilation to ob¬ 
servational gamma-ray photons. By using the current SGRBs data, the disk masses for the 
definite ranges of the BH parameters are shown in Section 3. Conclusions and discussion are 
in Section 4. 
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Model 


The neutrino annihilation luminosity L vP is a function of the BH mass Mbh, dimen- 


-- V"* / BH * “ --o "i -" 

mensionless viscositv parameter a , and mass accretion rate M (see, e.g., Popham et al. 

1999|; Rosswog et al. 

2003i Gue 

; al. 

2006|; Liu et al. 2007; 

Zalamea & Beloborodov 2011; 

Kawanaka et al. 2013 

Xue et al. 

2013 

; Leng & Giannios 201 

4). 


The analytical formula of L„ p is shown in m any previous works (e.g.. lFrver et al.lll999 


Zalamea fe Beloborodov 201 ll: Xue et al. 2013 ). Here we adopt the neutrino annihilation 
luminosity L uP given by Zalamea fe Beloborod ov (2 0111) . which is expressed as 


L,yi / 


5.7 x 10 
0 


52 - 4.8 

•''ms 


(M bh /M s ) 


- 3/2 


for M < M\ 


ign 


X 


(M/M® s" 1 ) 9 / 4 for M ign < M < M trap 
(M trap /M 0 s" 1 ) 9 / 4 for M > Mt rap 


erg s 


-l 


( 1 ) 


where x ms = r ms /r g , r ms is radius of the last (marginally stable) orbit, r g = 2 GMbh/c 2 is the 


rate if neutrino trapping events occur in the inner region of the NDAF (e.g.. 

Di 

Vlatteo et al. 

2002; 

Kohri et al. 

2005; 

Chen & Beloborodov 

2007; 

Liu et al. 

2012a; 

Xue et al. 

2013 

). Their 


numerical results depend on the viscosity parameter a and BH spin parameter a*. Addition- 


ally, the value of viscosity parameter a has little effects on L uD as long as M\ m < M < Mt. ran 

(Zalamea & Beloborodov 

2011 

), so a = 0.1 is adopted here. Furthermore, x ms can be ex- 

pressed as (e.g., 

Bardeen et al. 

1972; 

Kato et al. 

OO 

O 

O 

CN 

Hou et al. 

2014 

) 

*£ms — r[3 + Z 2 — \/(3 — 2fi)(3 + + 2 Z 2 )], 

(2) 


where 


Zx = 1 + (1 - a 2 ) 1/3 [(l + a*) 1//3 + (1 - a*) 1 / 3 ], 


(3) 


Z 2 = J3al + Z 2 . 


In comparison, Xue et ah ( 20131 ) also gave a similar analytical solution, i.e., L v 


(4) 

a M 2 - 17 , 


but the influence of the BH mass was not considered. 


Popham et al. (1999) and Liu et al. (120071 ) investigated the spatial distribution of neu¬ 


trino annihilation rate and found that nearly 60% of the total annihilation luminosity is 
ejected from the region r < 20 r g . In the studies on the vertical structure of NDAF model 
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(e.g., Liu et ah 2010, 2012a, 2013), we found that the half-opening angle of the disk is very 
large, > 80°, for the typical accretion rate, ~ 1 M 0 s -1 , thus the empty funnel along the 
rotation axis above the disk can naturally limit the opening angle of the neutrino annihilable 
ejection to produce the primary fireball. 

The fireball mean power outputting from the central engine if is a fraction of i.e., 

E = 7]L uD , (5) 


where ij is the conversion factor (e.g., lAlov et al.l 12005c iFan fe Weil 120 lit iLiu et al.ll2012bh . 
The output power can be written as 


• (1 + Z )(E-y iso + ifk,iso)^f 

E w- J 


271 


( 6 ) 


90 


where z is the redshift, E Jt i ao is the isotropic radiated energy in the prompt emission phase, 
£k,iso is the isotropic kinetic energy of the outflow powering long-lasting afterglow, T 90 can 
roughly be considered as the duration of the activity of the central engine, and dj is the 
opening angle of the ejecta. 


Hence, for the cases of M\ zn < M < M trap , we have the mean accretion rate (Fan fc We i 


2011 ) 


JiV ~ n in r *' 1 + (-^ 7 , iso,51 + -^k,iso,5l)^ 2 l4 / 9 _ 2 .1 

[ rfr 9 o, s J Xms m q 


M q s" 1 , 


(7) 


where £ k , iS o,5i = F; ki i SO /(10 51 ergs), A 7iis0i5 i = A 7iiso /(10 51 ergs), and T 9 o,s = T 90 /(l s). 
Furthermore, the disk mass is 

M Jilk « 0.12 | (E - , 'i'°’ 8 ‘ + - Ew ‘ ) ^ ]4 /9 (2BL)»/« ^ M 0 . (8) 

7] i + Z tkz 0 

According to the above equation, we can estimate the disk mass by using the obser¬ 
vational data. It should be noted that there exist some uncertainties, especially for the 
efficiency r] and the interval of the activity of central engine replaced by Tgo- There should 
exist an efficiency from the neutrino annihilation to the initial fireball, then to the jet ki¬ 
netic energy and r adiation, which is mainly related to t he energy, c ompo nents, and state 
of the fireball (e.g., Eichler et al. 1989 : Alov et al. 2005 b Alov et al. ( 2005 1 mentioned that 
the duration of the GRB event might be longer than the time interval of the activity of 
central engine if the radial expansion of the fireball is considered. In the fireball model, it 
is difficult to estimate the duration of such an expansion to the optically thin phase by the 
observational data unless the blaekbody component can be observed. It is conceivable that 
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the consequences of the use of r/ and Tgo would change the resulting disk mass to some extent 
although the exponents of rj and T 90 in Equation (8) are small. 

E liiso can be calculated by the observational data, which is written as 

£ 7 ,iso = 4ttD\F 1 / (1 + z), (9) 

where Dl is the luminosity distance, and F 7 is the fluence in the 15-150 keV for Swift events. 
Then D L is dehned as 

d l = { 1±^L i + ^) 3 + }- 1/2 dz’, (io) 

Jo 

here we employ a standard ACDM cosmology model with Om = 0.27, 12a = 0.73, and 
H 0 = 71 km s -1 Mpc -1 . Moreover, the mean isotropic gamma-ray luminosity is 

L'y, iso ~ 27 7i ; so (1 + z)/Tqq. (11) 




E\, is o and 9\ can be deduced from the modeling of the X-ray afterglow da ta. We take 
k)iso as ( Llovd-Ronning &: Zhang 2004 : Fan &: Piran 2006 : Zhang et al. 2007b ) 

1 + z, 


7 k,iso 


9.2 x 10 b2 RL\ /{p+2) ‘ 


)-l (2-p)/(p+2) 4(l-p)/(p+2) 


J A',46 V 2 ' B ~ 2 

x4 3p -2)/( p+2 )( 1 + y) 4 /(p+ 2 ) ergs, (12) 

where R ~ (tn/T 90 s ) 17ee L 6 is a factor that accounts for t he energy loss during the decelera¬ 
tion f ollowing the prompt gamma-ray emission phase (e.g.. lSarilll997l:lhlovd-Ronning fc Zhang 


20 04), e e ,o.i = e e /0.1 is the fractions of shock energy given to the electrons, € 3 , 0.1 = Cb/0.01 


is the fraction of energy in the magnetic field, tn = 7/(11 hours) and t, ( \ = 7/(1 day) are 
the time of observation, Y is Compton parameter, p is the energy distribution index of the 
shock-accelerat ed electrons and c a n be fitted by t he observed photon index in the X-ray 
spectrum (e.g., IZhang et al.l 120061 : lOao et al.l 120131 ). and Lx ,46 = L x /{ 10 46 erg s x ) is the 
isotropic X-ray afterglow luminosity. Here we take the X-ray luminosity at 11 hours since 
the burst triggers, which can be written as 


L x — 4nD L F x , 

where F x is the X-ray flux of the afterglow recorded by satellites. 


(13) 


Furthermore, the relation between the opening angl e and the jet break time is given by 


(e-g., 


Sari et al.lIl999l: iFrail et al.ll200ll: iFong et al.l 120121 ) 


0j ~ 0.076 


L \3/8/ 1 F Z \-3/8f 


n 


_ R / s??- 1 / 8 

,— 3 / ^k, iso, 51’ 


(14) 


1 day' v 2 ' v 0.01 cm -3 ' 

where t- } is the jet break time in the X-ray afterglow phase of GRBs, and n is the number 
density of the burst circumstance. 










































7 


3. Results 

The relations among the observational data of the prompt emission and afterglow in 
SGRBs, disk mass, and BH parameters are established by Equations (8-14). If the reasonable 
ranges of the BH parameters are given, the limits of the disk masses corresponding to the 
certain SGRBs can be resolved. 


3.1. Data of SGRBs 


Berger (2014) mainly reviewed the progresses of SGRBs in the theories and observations, 
including the afterglow and host galaxy observations, the properties of the circumburst 
environments and their progenitors. There are 70 SGRBs with a substantial fraction of 


afterg low detections in the eight-year period from January 2005 to January 2013 (IBerger 
2014), with the addition of G RB 130603B (e.g ., Berger et al. 2013 : Tanvir et ah 2013), which 
is associated with a kilonova (lLi fe Paczvnski 1998]). As shown in Table 1 of Berger) (2014), 
there are 27 SGRBs with the authentic X-ray detections and known redshifts discovered by 
Swift satellite except for GRB 050709 by HETE-2. 

Moreover, we find 4 SGRBs with the X-ray detections and known redshifts triggered 
after GRB 130603B, i.e., GRBs 131001A , 14 0622A, 140903A, and 141212A, whose data are 
from the UK Swift Science Data Centre ( Evans et al. 2009 ). So totally 31 SGRBs are listed 
in Table 1. F or each SGRBs, we fit the photon index with the data of X-ray spectrum to 


deduce p (e.g. JZhang et al.ll2006l: iGao et al.ll20131) . Their durations T 90 , redshifts z, gamma- 


ray fluences F 7 and X-ray fluxes at 11 hours since trigger Fx (11 hours), and the observed 
spectral index f3 are displayed. If we take e e ~ 0.1, 6 b ~ 0.01, Y ~ 0, rj — 0.3, and fitted p 
and given tj, then E 7;iso , E^ iso , and 6- } can be solved, and we further obtain the ranges of the 
disk masses. 

It is worth noting that the most difficult problem is the estimation of the jet opening 
angle 9 } because of the faint and restricted observations of SGRB afterglows. So far there 
are three scenarios, i.e ., (1) a fe w credible detections of jet brea k, such as in G RBs 051221A 
( Soderberg et al.ll2006 ). 090426 ( Nicuesa Guelbenzu et al. 2011 ). and 130603B ( Fong &: Berger 
201 3lf : 12) several meaningful lower limits on jet op ening angles, s uch as in GRBs 05 0724 
( Grupe et al. 2006 ). 111117A ( Margutti et al. 2012 ). and 120804A ( Berger et al. 2013 ): (3) 
no break in X-ray lightcurves of some SGRBs. We cite the data of the jet opening angles 
or their lower limits in the above references for former two cases as shown in Table 1. For 


the third scenario, we set the lower limit of 9j > 0.05 (IFong et al. 2012). For the last 4 


SGRBs we collected, the jet break time, ~ S.SOj^so ks, is found in GRB 140903, then we 



























































can estimate the opening angle by Equation (fl4l) with n = 0.01 cm 3 , and the other three 
are set by the lower limit as discussed above. 


3.2. Disk masses of SGRBs 


Figure 1 shows the ranges of disk masses Mfcsk of the different SGRBs with the isotropic 
gamma-ray luminosity L 7 j i SO for varying BH mass Mbh from 2.7 M 0 to 10 M 0 and BH spin 
a* from 0 to 0.99. Three vertical lines correspond to = 0.2 M 0 , 0.5 M 0 , and 1 M 0 , 
respectively. It is seen that M disk has a wide distribution, from about 6 x 10~ 4 M 0 to about 
7.6 Mq, as well as shown in Table 1. The accretion rates corresponding to the minimal disk 
masses are checked, which are in the suitable ranges, i.e., M; gn < M < M trap . There is no 
statistical correlation between the disk mass and gamma-ray isotropic luminosity, because 
the energy coming from the accretion powers all radiative processes of GRBs, mainly in 
gamma-ray and X-ray bands. The energies of the X-ray afterglows are frequently larger 
than those of prompt emission as displayed in Table 1. Obviously, the disk mass is primarily 
calculated by the output energy of GRBs, opening angle of the jet, and BH characteristics 
as shown in Equation (8). There exists a difference of several orders of magnitude between 
the minimal and maximal disk mass, which means that the BH characteristics are the major 
factors on the disk mass. 


As shown in the figure, the maximal disk mass of GRBs 050724, 051221A, 070714B, 
070809, 090426, 111117A, 120804A, and 131001A are larger than 1 M 0 , and that of most 
other SGRBs in our sample are larger than 0.2 Mq, which indicates that the extreme BH 
spin parameters and small BH mass are required. In binary NS merger events, the BH mass 
is naturally less than the total mass of the binary, i.e., ~ 4 Mq, which is described about 
~ 2.7 Mq in further simulations. In BH-NS binaries, the BH is origin from its progenitor 
star, and its mass should also be a stellar-mass order, e.g., ~ 10 Mq. Moreover, _the_BH 


spin p aram eters are also related to their progenitors. In some discussions (e.g., iLee et ah 
2000b; Ruffert & Janka 2001), the rapidly rotating BHs (a* > 0.5) are inclined to exist in 


the SGRB centers in contrast to the BHs in the LGRB centers. 


In order to embody the effects of the BH characteristics on the disk mass, Figure 2 dis¬ 
plays the distributions of the disk masses M &isk f° r the different typical BH masses and spins, 
which are set to (M B n/M Q , a*) = (3, 0.5), (3, 0.9), (10, 0.5), and (10, 0.9), corresponding 
to (a-d), respectively. It is easy to find that the disk mass of most SGRBs are safe below 
0.2 — 0.4 Mq, and sporadic cases are beyond the limits, especially in the case of Figure 2(c). 
Additionally, by comparing these four cases, we notice that the spin parameters are more 
effective than the BH masses on the values of the disk masses. Even for the case of Figure 
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2(b), there still exists one SGRB, whose disk mass is larger than 0.45 M & . Those massive 
disks, > 1 M 0 , may exist in the centers of collapsars, which are considered as the origin of 
LGRBs. Actually, Lazzati et al. (2010) proposed that the off-axis jets from collapsars could 
power SGRBs. 

Three factors remind us that the results of the disk masses are the lower limits at the 
most, which are as follows. (1) We have to calculate the disk mass using the lower limit of 
in most SGRBs as shown in Table 1. It is easily conceivable that the real requirements 
of the disk masses are much larger than the the present results if the precise value of 
is considered. (2) Some powerful SGRBs with unknown redshift, such as GRBs 060121, 
060313, and 111121A, shown in Table 1 of Berger (2014), may require more massive disk 
than SGRBs in our sample if they are also origin from the BH hyperaccretion systems. (3) 
The powerful X-ray flares have been extensively observed in the afterglow phase of GRBs, 


which are considered t o originate from the re-ignition of the central engine (e.g., iLiu et al. 


2008 


Luo et al.l 20131 : Hou et al. 2014). This means that the remanent matters from the 
massive disk are needed to maintain the explosion of X-ray flares. However, we use Tqq to 
replace the duration of the activity of the central engine, which may generally enlarge the 
disk mass in the calculations. Although these influences and some uncertainties may exist, 
we consider that our results can still reflect the deficiency of the neutrino annihilation process 
to power SGRBs. 


4. Conclusions and discussion 


The progenitor of SGRBs is considered to be a compact binaries merger event. After 
merger, a stellar-mass BH surrounded by an NDAF will be formed in the central SGRB and 
neutrino-antineutrino annihilation above the disk may power SGRBs. The total mass of two 
compact stars limits the mass of the system consisting of a BH and an NDAF. In this paper, 
we focus on a question, i.e., can annihilations of neutrinos from NDAFs owning such masses 
power all the observed SGRBs? The calculations show that the disk mass of a certain SGRB 
mainly depend on its output energy, jet opening angle, and central BH characteristics. Even 
for the extreme BH parameters, there still exist some SGRBs requiring the massive disks, 
which approach or exceed the limits in simulations. 


Besides magnetar model, for BH hyper accretion system, w e suggest t hat there m ay ex ist 
an alte rn ative magnetic or i gin of SGRBs, i.e., BZ proces s (e.g., iPopham et al.lll999l;iLee et al. 


2000 al]bl: iDi Matteo et al.l 120021: iKawanaka et al.l 120131) or episodic magn e tic re connection 


pr 

J i 


(Yuan fe Zhangi 2012), to replace neutrino annihilation. Kawanaka et al. (120131 ) presented 
that the luminosity powered by Poynting-dominated jet is more qualified for the requirement 
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of GRBs than neutrino pair annihilation. Yuan & Zhang (2012) investigated that the closed 
magnetic field lines continuously emerge out of the accretion flow. Since the shear and 
turbulent motion of the accretion flow, the line may form the flux rope. When a threshold is 
reached, the system loses its equilibrium and the flux rope is thrust outward, then an episodic 
jet occurs. This mechanism can also power enormous energy to trigger GRBs. In addition, 
if these magnetic origins really exist in the center of GRBs, the polarization effect should be 
observed in the prompt emission or afterglow of GRBs. Actually, the linear polarization in 
the afterglow of LGRB GRB 120308 has been detected I Mnndell et al. 2013 b which indicates 
that large-scale magnetic fields may be dominant in the GRB jets. But now we do not know 
whether the same situation exists in SGRBs. 

Otherwise, there are some mechanisms, such as magnetic coupling from the BH horizon 
to the inner region of the disk (Li 2000), can effectively transfer the angular momentum 
and rotational energy of the BH to heat the inner region of the disk, then l arger num bers 


of neutrin os r adiate from the disk to produce the primordial fireball (e.g., [Lei et al. 2009 


Luo et al.ll2013lh Besides, the vertical advection (or convection) i s considered to widely exist 
in the slim disks and NDAFs (jji an g et ah 20.14; [Liu e t a l. 2015,), which is another possible 
mechanism to increase the neutrino emission rate. The scenario for NDAFs follows, the ver¬ 
tical advection (or convection) caused by magnetic buoyancy can much effectively transport 
energy to the disk surface, and also suppress the radial advection, thus the neutrino lumi¬ 
nosity and annihilation luminosity are dramatically increased. This mechanism is conducive 
to achieve the energy requirement of GRBs. 
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Table 1: Data of SGRBs 


GRB 

T 90 

00 

z 

F 1 

(10 7 erg cm 2 ) 

Ex (11 hours) 

( 10 —14 erg cm -2 s — 1 ) 

Photon index 

^ , y,iso 
( 10 5 ergs) 

^k.iso a 
( 10 5 ergs) 


(rad) 

M disk 6 

(M^) 

050509B 

0.04 

0.225 

0.23 

< 1.95 

l g+U.5 

0.0027 

0.055 

> 

0.05 

0.0006-0.028 

050709 

0.07 

0.161 

4.0 

1.92 

~ 2 

0.023 

0.016 

> 

0.26 ( 1 ) 

0.003-0.14 

050724 

3 

0.257 

6.3 

9.55 

i-ssioR 

0.1 

0.27 

> 

0.35 (2) 

0.084-3.93 

051210 

1.3 

1.3 

0.83 

< 2.7 

r\ iyoTO-48 

° —0.41 

0.36 

2.38 

> 

0.05 

0.016-0.76 

051221A 

1.4 

0.5465 

12 

108 

o 0 Q + 1 '*'•*' 9 

0.09 

0.92 

12.6 


0.12 (3) 

0.093-4.37 

060502B 

0.09 

0.287 

0.4 

< 1.47 

9 -1 e+i-U < 
2 - iD -0.58 

0.012 

0.12 

> 

0.05 

0.0013-0.062 

060801 

0.5 

1.130 

0.81 

< 0.98 


0.27 

0.71 

> 

0.05 

0.0063-0.30 

061006 

0.4 

0.438 

14 

22.7 

i-sefo.lS 

0.67 

3.14 

> 

0.05 

0.013-0.59 

061201 

0.8 

0.111 

3.3 

19.2 

1 54 + 0 - 17 

0.17 

0.01 

0.07 


0.017 (4) 

0.0014-0.067 

061210 

0.2 

0.409 

3.0 

13.6 

2 60”*” 1-92 
^.OU_o ?! 

0.12 

0.86 

> 

0.05 

0.0048-0.22 

070429B 

0.5 

0.902 

0.63 

11.3 

o fiq+bl 8 
^.oy_o 56 

0.13 

4.51 

> 

0.05 

0.013-0.63 

070714B 

2.0 

0.923 

7.2 

6.30 

1 96+ U-1 ^ 

1,yD —0.15 

1.61 

2.32 

> 

0.05 

0.027-1.25 

070724A 

0.4 

0.457 

0.30 

12.8 

1 46"*" 0 ' Sti 
—0.25 

0.016 

0.99 

> 

0.05 

0.007-0.33 

070729 

0.9 

0.8 

1.0 

< 4.71 

i c+ 0.6 

10 -0.3 

0.17 

1.32 

> 

0.05 

0.012-0.54 

070809 

1.3 

0.473 

1.0 

53.0 

1 ^Q+0.14 

i.3y _ 0 12 

0.056 

3.91 

> 

0.05 

0.024-1.15 

071227 

1.8 

0.381 

2.2 

3.20 

o i q+0.41 
2 ' iy -0.35 

0.08 

0.25 

> 

0.05 

0.01-0.47 

080905A 

1.0 

0.122 

1.4 

< 6.7 

1 54 + 0 - 22 
a -°^_0.14 

0.005 

0.024 

> 

0.05 

0.0027-0.13 

090426 

1.2 

2.609 

1.8 

26.3 

2 nQ +0 - 1 - 6 
Z.U3_o 15 

2.84 

135 


0.07 (5) 

0.093-4.35 

090510 

0.3 

0.903 

3.4 

5.04 

1 70”*” 0i2 

A - ' U - 0.12 

0.73 

3.07 


0.017 (4) 

0.0035-0.17 

090515 

0.04 

0.403 

0.21 

< 8.43 

9 7 oT 1-20 
<O-0.?7 

0.008 

0.62 

> 

0.05 

0.0016-0.075 

100117A 

0.3 

0.915 

0.93 

< 2.50 

2 74 +°- 36 

1 ^—0.31 

0.20 

1.10 

> 

0.05 

0.0057-0.27 

100206A 

0.1 

0.408 

1.4 

< 1.07 

o O "*” 0,8 
zu -0.7 

0.058 

0.073 

> 

0.05 

0.0013-0.062 

100625A 

0.3 

0.453 

2.3 

0.395 

9 q + 0.5 

0.12 

0.093 

> 

0.05 

0.003-0.14 

101219A 

0.6 

0.718 

4.6 

2.00 

1 44+0-27 

1 -^-0.25 

0.62 

0.45 

> 

0.05 

0.008-0.38 

111117A 

0.5 

1.3 

1.4 

3.21 

2 in+0-39 
z,iU -0.32 

0.62 

3.77 

0 . 

105 ( 6 ) 

0.023-1.06 

120804A 

0.81 

1.3 

8.8 

58.6 

2 in+ 0-22 

^• iU -0.14 

3.88 

56.9 

> 

0.19 (7) 

0.16-7.59 

130603B 

0.18 

0.356 

6.3 

60.0 

2 00”*" 014 
^.UU _ 0 13 

0.20 

2.80 


0.07 ( 8 ) 

0.01-0.48 

131001A 

1.54 

0.717 

2.8 

14.7 

1 9l + 018 
x ‘ y ‘ L —0.18 

0.37 

5.41 

> 

0.05 

0.029-1.37 

140622A 

0.13 

0.959 

0.27 

17.0 

1 cc+0.67 

i.oo_o 28 

0.065 

9.77 

> 

0.05 

0.0087-0.41 

140903A 

0.30 

0.351 

1.4 

124.7 

1 rq+ 0.22 

x.o »_0 20 

0.043 

6.15 

0 . 

023 (9) 

0.0067-0.31 

141212A 

0.30 

0.596 

0.72 

2.50 

2 o + 0,8 
^• U -0.5 

0.066 

0.38 

> 

0.05 

0.0039-0.18 


Notes: 


a The parameters are calculated by Equation J 12 1 with e e 0.1, eg 0.01, and Y ~ 0. 

^ The ranges of are estimated by Equation f 8 l with rj = 0.3, varying M-q from 2.7 Mq to 10 Mq, and a* from 0 to 0.99. 

References: 
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2006J; (3) Soderberg et al. 20061; (4) De Pasquale et al. 20ld; iNicuesa Guelbenzu et al. 20121; (5) 

Nicuesa Guelbenzu et aL||201l]; ( 6 ) 

Margutti et al. 20121; (7) iBerger et al. 20^1; ( 8 ) |Fong & Berger 20^; (9) The opening angle of GRB 


140903A is determined by Equation 1141 with the data from the UK Swift Science Data Centre {Evans ct al 112009ft and n ~ 0.01 cm 8 . 



















16 



Fig. 1.— Ranges of disk masses Mh j H y of different SGRBs with isotropic gamma-ray lumi¬ 
nosity L 7) i S o for varying BH mass Mbh from 2.7 M 0 to 10 M 0 and BH spin a* from 0 to 
0.99. Three vertical lines correspond to Mdi S k = 0.2 M 0 , 0.5 M 0 , and 1 M 0 , respectively. 
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Fig. 2.— Distributions of the disk masses M c 



for different typical BH masses and spins 











